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Motivated by phenomenological questions in quantum gravity, we consider the propagation of a 
scalar field on a random lattice. We describe a procedure to calculate the dispersion relation for the 
field by taking a limit of a periodic lattice. We use this to calculate the lowest order coefficients of 
the dispersion relation for a specific one-dimensional model. 
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I. INTRODUCTION 

7-H ' 

In the present article we will explore the propagation of a scalar degree of freedom on a one-dimensional lattice. The 
lattice is constant in time, but the geometric information about the lattice - encoded in a metric-like function - is 
I allowed to vary spatially. [13] In fact we are most interested in the case where the geometry fluctuates in a stochastic 
way along the lattice. We call this a random lattice. The dynamics we consider for the scalar field is that of a massless 
Klein-Gordon field, discretized on a lattice with variable lattice spacing. Thus, in the limit of the lattice spacings 

■ going to zero, one expects to obtain the well known continuum model. Here we are however interested in the model 
t-H ' with finite lattice spacing, but for excitations of the field that are on very large scales as compared to the average 

lattice spacing. In this regime, the system should show propagating waves, governed by a dispersion relation. The 
' question we pose to ourselves is how this dispersion relation depends on the microscopic parameters of the theory. 
Q" 1 This is very much in analogy to the question of how to determine the speed of sound, and other macroscopic properties 
of a (possibly amorphous) solid, from the microscopic properties of its constituents. We are, however, not motivated 
by solid state physics, but by quantum gravity. 
I— ~~ 1 ■ The incorporation of the quantum nature of gravity is likely to change the structure of space-time on small length- 
| scales. Different such structures have been investigated, some motivated by full-fledged attempts to quantize gravity, 
some as toy models or as effective descriptions. One paradigm is that spacetime is fundamentally discrete, such as 
^) in the causal set approach [lj] , or in semiclassical considerations of loop quantum gravity . Another paradigm, 
pioneered by Wheeler, has spacetime smooth, save for localized defects 5]. Our model here is an instance of the 
former, but could maybe also be read as a model of space-time defects in the limit where their density becomes very 
high. The model is very similar to the ones considered in loop quantum gravity 043] , and references therein (and 
some intermediate results of the present work were already contained in |6|). There, dispersion relations were also 
t-H ' calculated, but the calculations involved simplifications, effectively replacing a random lattice by a regular one, the 
parameters of which were obtained by averaging. In the present article we go beyond this, by obtaining exact results. 
We should point out that because we are working with a free massless field in two dimensions, the corresponding 
continuum model would be symmetric under the exchange of space and time coordinates. Therefore, our model 
is mathematically equivalent to one in which spatial geometry is smooth and homogeneous, but time is a discrete 
variable. For simplicity, we will work in the "discrete space, continuous time" picture, but we will also give the results 

■ for the opposite case in an appendix. 

Since we are working with random lattices, our results are stochastic in nature. So while we said that we would like 
to compute the coefficients in the dispersion relation, what we actually compute are expectation values (and in one 
case the variance) of such coefficients. 

Fields on random lattices have become a valuable tool in lattice gauge theory, starting from the pioneering works 
0,9- We are, however, not aware of explicit results on the dispersion relations of these fields at finite lattice spacing. 
This may have to to with the fact that in this context, the random lattices are a tool to obtain statements on the 
continuum limit of the theories. 

Effects of coupling to random fields on the dispersion relation have already been studied in detail in Q and [IcJ • In 
these works there is no discreteness of space-time. Rather, the random fields are an effective description of the effects 
of the CPT anomaly, caused by a nontrivial space-time topology on small length scales. While these models are vastly 
more physical, and can hence be used to obtain bounds on small-scale space-time structure, one can nevertheless 
compare them to our model on a mathematical level. The results bear some intriguing similarities but also some 
differences. These are discussed in a bit more detail at the end of this article. 

It should be said that our model is rather unphysical, in that it is 1+1 dimensional, and in that the field content is not 
realistic at all. What we hope to have accomplished, is to show how a calculation of phenomenologically interesting 
data from a model involving discrete, randomly fluctuating space-time can be accomplished. Our calculations here 
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can almost certainly be generalized to more realistic models (using, for example, the Voronoi-construction to 
obtain a random lattice in arbitrary dimension). They would merely be more cumbersome. 

The structure of the article is as follows: In the next section we will specify the details of the model. Section IIIII 
explains the calculation of the dispersion relation. We end with a discussion of further prospects in section IIV1 Two 
appendices contain some of the longer calculations, and the last appendix contains the dispersion relation for the case 
of continuous space and discrete time. 

II. THE MODEL 

We consider a bosonic field propagating on discrete space and continuous time. (The case of discrete time and 
continuous space is briefly discussed in appendix O) To keep things simple, we work in 1+1 dimensions. The field is 
a function on the space IxZ. It will be denoted by <j) n {t). The geometry of space is encoded in a time- independent, 
positive function g n . Often it is also useful to re-express g n as g n = l~ 2 . At this point we will make no further 
assumptions on g n , but we will later assume that it is random in a certain specific sense. 
The action for our model reads 

S = f dtJ2l[<P 2 n-9n(d + <P) 2 n ]. (1) 

Here d + (f> is the forward discrete derivative, (d + (f>) n = cj> n +i — 4>n- It w &s used for simplicity. A more symmetric 
derivative could be used in its place. We will also use its adjoint d~ with respect to the sum over Z. S can be viewed 
as a discretization of the action for a free massless scalar in the continuum, but ultimately an action of this form may 
be derived from a theory of quantum gravity, and thus regarded as fundamental. 
The equations of motion are 

4> n + A g (j) n = 4> n - g n <p n +l - 9n-l4>n-l + (ffn-l + ffn+l)0n = (2) 

where we have introduced the discrete, positive definite Laplacian 

A 9 = d-g n d+. (3) 

For completeness, we also state the Hamiltonian 

H =\H<+9n{d + ^)l- (4) 

n 

It is of the form that is obtained under certain assumptions in loop quantum gravity [H, Q • 

The form of the function g n is obviously vital for the definition of the model. For our consideration, three cases are 
of interest: 

1. g n is constant 

2. g n is periodic, with a certain period TV in Z. 

3. g n is obtained as an instance of a random process. 

Case 3 is the one we want to consider. The simplest situation would be that the value of g ni would, for any given n, be 
determined by the sampling of a certain random variable g. Another way to state this is that there are independent, 
equally distributed random variables g„ and the function g n is a sampling of these. More complicated situations (for 
example correlations between the g n ) are also conceivable. 

As we will see, it will be important to have information about the distribution of values {g n } in the sampling of {gn}. 
Of particular importance for us will be the mean of the l n , as well as some related quantities. Such information may 
in principle be obtained through the law of large numbers, or a central limit theorem. 

To be concrete let us fully specify one simple model: To make the description simple, we will not specify the distribution 
of the variables g„, but those of the 1„. We denote expectation values by angular brackets (•). Let all the random 
variables 1„ be independent Gaussian distributions with first and second moments 

(1„)=Z, ((l n -l) 2 ) = d 2 . (5) 

Returning to the list of cases from above, the first case is that of an equidistant lattice. It is immediately solvable. 
The second case is still solvable and we will use it in order to analyze the case we are really interested in. We will 
see that we can obtain the lowest order terms of the dispersion relation for case 3 in a limit of case 2. This will be 
explained in the following section. 
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FIG. 1: Optical and acoustic branches in the dispersion relation 



III. DISPERSION RELATIONS 



The problem of wave propagation on a random lattice is certainly a complicated one. In fact, we will see indications 
that it is not even always well defined, i.e. that there is not always a long wavelength limit in which something 
resembling plane waves propagates on the lattice. The question is how to identify the cases in which the problem is 
well defined, and how to extract characteristic long wavelength quantities. We will not address these problems in all 
generality. But we will obtain a formula that, in a certain limit, gives the first few terms in the dispersion relation for 
the field. This limit is by no means well defined for all random lattices. If it is ill-defined, this is a strong hint that 
there is no regime in which the lattice supports propagating waves. 

To start our quantitative discussion, we consider the case of a regular lattice. If all the U are equal (to /, say), it is 
easy to solve the equations of motion of ([2]). The solutions are "plane waves" 

n (* ) fc) = e < ( Wn —W*) ) u J 2 (k) = ^(l-cos(kl)) = k 2 -^k 4 + 0(k 6 ). (6) 

For the general case on the other hand, with generically all Zj different, it is not possible to explicitly write down 
any solution to the equations of motion. The analysis we are aiming at in the present section lies somewhere in- 
between these two extreme cases. We will make an assumption on the k under which we are able to treat the system 
analytically and try to remove it at the end of the analysis: Let us assume that the system is periodic with N £ N the 

(z) 

length of period. More precisely we assume that g n +N = 9n for all n G Z. We introduce the notation <p„ = <j) n +zN 
with n £ {0, 1, . . . , N — 1} and make the Ansatz 

N-l 

(t) = Cn expi(Lzk -tot), L=J2l n . (7) 

n=0 

This Ansatz turns the equations of motion © into an eigenvalue problem for c and w: (0 is a solution iff 



Mc = uj 2 c where M = 



(9N-1+90 -go ... -gw-ie \ 

-.9o 9o + 9i -.9i ... 

-.91 5i + .92 -52 

V-5JV-ie~ lfcL ... -5jv-2 9N-2+ 9N-J 



(8) 



The eigenvalues loq . . .u>n-i represent the different branches of the dispersion relation. We will presently see that 
there is one branch, denoted uj ac in the following, with w ac (fc) — > for k — > 0. Following the custom of condensed 
matter physics, we call this branch acoustic in contrast to the optical branches nonzero at k = 0.[l4j The situation is 
sketched in figure Q] As we are interested in the low energy (i.e. small lu) behavior of the field, the acoustic branch 
is the relevant one for our purpose and we will compute its small k behavior in the following. Let us start by making 
the Ansatz 

Lul c (k) = w(l) \k\ + w(2)k 2 + w(3) \k\ 3 + . . . , (9) 
explicitly forcing w ac (0) to be zero. Accordingly, we expand det(M — w 2 I): 

N-l oo 

det(M-^ 2 I)= u 2l Y,w{i,j)\k\ 3 . (10) 

i=0 3=0 
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We want to determine the coefficients w(i) of ©. The coefficients w(i,j) of the expansion (TTOf on the other hand 
can be explicitly calculated. The calculation is cumbersome, so we have relegated it to appendix [X] We find 
Proposition III.l. For M_ of the form (jHJ), 

det(M-^ 2 I) = - 2g ... g N ^i(l - cos Lk) 

AT-l 

2 Ng . ..g N -i Si* 



N-l 

(11) 



+ W 4 9o . . . 5JV-1 51 0' - i N - U - *)]9i l 9'j 1 

0<i<i<JV-l 

+ 0(oj 6 ). 

Froom this, the coefficients w(i,j) can be read off. By solving the eigenvalue equation det(M — cj 2 I) order by order, 
we then obtain the low order coefficients of the dispersion relation ([9]) . We use the shorthands cy = (j — i)[N— (j — i)] 
and 



1 N - 1 

f = Jf E /« (12) 



i=0 

for the average of some quantity over the period of the lattice. Then we have 



u 2 ( k ) _ L% l 7,2 , f 7 - : ^- - r -<- y - "j ' 1 Li \ k 4, (k 6 ) 

u, ac (fc) - N2 _fc + I 6 _ I A + O (A j 



This is a remarkable formula, and one of the main results of the present work. It is an exact result for the lowest 
orders of the dispersion relation of the field 4> propagating on a periodic lattice. As such, it reproduces the elementary 
result ((6]) for the case of the regular lattice, by setting l n — I for all n and using 



E 

0<i<j<N-l 

We note that the square of the velocity of propagation, 



^ = ^NHN 2 ~ 1). (14) 



t 

c 2 := L (15) 
P 

is neatly expressed in terms of averages over lattice spacings. The formula for the next order coefficient, 

1 f T 2 1 2 
^^J-^V^-iL (16) 

\ I i<j 

is certainly more complicated, but also given in terms of such averages. This is more than just aesthetically pleasing: 
It suggests that (fT3")) may survive the large- TV limit. Moreover, that when interpreting the U as given by a sampling 
of some random process, such averages may be expressible by averages under this process, through some ergodic-type 
results. 

Let us come back to our original motivation: We were interested in the case of a periodic lattice, because, for the 
l n determined by some random process, the limit N — > oo, with the average length (lj) = / held fixed, gives an 
infinite random lattice. Thus we would like to consider the dispersion relation (|13p in that limit. We see, however, no 
possibility to discuss this limit without fixing a distribution for the U. Central limit theorems would make statements 
about the expectation value of the mean 1 and functions of it in the limit N —> oo for large classes of distributions. 
They can similarly be used to obtain statements about the mean of the squares, l 2 , in this limit. But in our case, we 
need to investigate a function of both, I and P in the case of the first term in the expansion (1131) . The second term in 
(1131) additionally depends on N and on a curious linear combination of l 2 l 2 , thus further complicating the situation. 
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FIG. 2: Theory and numerical simulation for the coefficient c 2 : For I = 1, we plot the theoretical result (|18[) as a function of 
d (continuous curve) together with two numerical simulations (N = 100: small dots, N = 1000: large dots). Convergence is 
nicely visible. 



Let us therefore consider the situation described in ([5]) in which the k are independently Gaussian distributed, with 
average {h} — I and variance 

if _ = d 2 . What we want to do is to compute expectation values of the coefficients 
(|13[) under this distribution, in the limit N —> oc. 

Let us first consider the coefficient c 2 of the |fc| 2 -term in (|13[) . In appendix [Bl we show that while we can not compute 
its expectation value directly, we can expand it in powers of d/l, compute the expectation values of the first few terms, 
and take N to infinity. The result 

1 d2 dA d6 

- 1 -¥ + ¥-¥ + -- (17) 




strongly suggests, that 




l 2 



(18) 



For independent confirmation we have also checked this result on a computer, by calculating instances of c 2 for 
randomly generated sets {h}. Formula (|18|) describes the numerical results very well. An example is given in figureO 
Furthermore, we are able to calculate the variance of c 2 , again as a series in powers of d/l, and take the N — > oo limit 
term-wise. It turns out that the variance vanishes in this limit: 

lim (c 4 ) - lim (c 2 ) 2 = 0. (19) 

Details are again given in appendix |5] 

The situation is more complicated with the term proportional to |fc| 4 . Expansion in d/l (see appendix) gives 



which does not converge term by term in the limit N — > oo. But this does not say that the limit does not exist. In 
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FIG. 3: Theory and numerical simulation for the coefficient I 2 : For I = 1, we plot the theoretical result (|21[) as a function of 
d (continuous curve) together with two numerical simulations (N = 100: small dots, N = 1000: large dots). Correspondence 
between theory and numerics is visible for very small values of d only, and convergence apparently gets worse for larger N . 

fact, if one would extrapolate from the result for the coefficient c 2 , 

lim (P) = lim 

AT ^ \ / JV-^~ 



JV- 




1<J 



i 2 



(21) 



1 



I 4 



12 I 2 + d 2 
= -U 2 1 



(N 2 



12 I 2 + d 2 



N 



12 1+- 

one expects the result to be finite. Numerical results for small d/l indeed point towards convergence (figure [3]), but 
for larger d/l, convergence in the limit N —> oo could not be established. It may be there but too slow to be seen. 
For completeness, we have also considered the case of discrete time and continuous space in appendix[Cj The resulting 
dispersion relation is 



1 f 2 

uj 2 {k) = -U 2 - -T^k 4 
w c 2 c 6 



(22) 



However, we have not considered the expectation values of the coefficients in any detail. 



IV. DISCUSSION AND OUTLOOK 



In the present paper we have calculated the (expectation values of the) first coefficients in the dispersion relation 

Lo 2 (k)=c 2 \k\ 2 +f\k\ 4 + ... (23) 

for a scalar field propagating on what one could call a Gaussian random lattice in 1+1 dimensions. The first term, 
also setting the phase velocity, is unitless in our convention, thus the natural value is 1. We have seen 



lim (c 2 ) = — 



1 



P 



(24) 



i.e. it is in fact very close to 1, as long as the variance of the lattice spacing is small. We have also presented very 
strong evidence that its variance is strictly zero for an infinitely extended lattice, thus one does not even have to 
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take into account that in principle it has a probabilistic nature. In more realistic models, this coefficient would be 
very interesting for phenomenology: As was demonstrated in experiments and observations are very sensitive to 
differences of propagation speed of the fields in nature. The second term had a more complicated structure, and 
convergence in the desired limit could not be demonstrated to complete satisfaction. We have conjectured the limit 

lim (i 2 ) = -LJ?. (25) 

Its scale is set by the average lattice spacing squared, I 2 , as long as the variance of the lattice spacing is small. We 
have also briefly discussed a model in which the roles of space and time are, switched, i.e., time is discrete and space 
is continuous. There we see potential problems with causality, with phase and group velocity increasing due to the 
lattice effects. 

An interesting observation was made by a referee of a draft version of this article: The first two terms of the dispersion 
relation for the field on a regular lattice © can be turned into the first two terms of the dispersion relation on the 
random lattice, by rescaling k — >■ kl and dividing by (I 2 ) in ([5]). It is very well possible, that this procedure gives the 
correct dispersion relation on the random lattice to all orders. This does not detract from the value of the results 
presented in the article: While there may be a shortcut, it is still neccesary to show that the shortcut indeed provides 
the right results. In fact, one big motivation for the article was that several different "shortcuts" were used in the 
literature on quantum gravity phenomenology, with different results. We should also point out that, despite the fact 
that the dispersion relation of the random lattice seems to be obtainable from that of a regular lattice by the above 
scaling operations, it is not that of such a regular lattice. 

Similar, but much more detailed, calculations were carried out in 0, [l^ for random fields coupled to scalar and 
electromagnetic fields in physical space-time dimensions. Since these were continuum models, one expects differences 
in the results: By reformulation our model in terms of a continuum field, one can show that in the equations of 
motion, there are correction terms of arbitrary derivative order, as compared to the equations of motion for the free 
scalar field. In contrast, the fields of [9, [l(| receive only corrections in terms of first derivatives of the field. Still the 
results on the dispersion relation are similar: The phase velocity receives a correction downward, and the coefficient 
of fc 4 is negative, thus insuring causal propagation. Moreover, the mathematical structure is similar: In both cases, 
the second moment of the random field enters the phase velocity, and the coefficient of fc 4 is given by moments the 
autocorrelation of the random field, the first moment in case of 0, the second moment (of 1%) in the present 
model. 

Besides finishing the discussion begun in the present work, by showing convergence or divergence of the second 
coefficient t 2 , and studying the complete model for more generic distributions for the random lattice, it would be 
interesting to study the nature of the eigenvectors c, at least in low order in k. For example: Does the one for the 
acoustic branch really look like a plane wave, at least for large Nl Obvious further goals are the extension of the 
formalism to vector and fermionic fields and to physical dimension. At that point, one would be in a position to 
compare the models with experiment and put bounds on parameters of, for example, loop quantum gravity. 
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Appendix A: Proof of proposition 1111,11 



It is an elementary combinatorial fact that 

JV-l 



t-^ , / sum over all (N — i) x (N — i) \ 

det (M-u 2 I )= V (-1 ) 1 uj 2i \ . . , , , r t r \. Al 

— \ principal sub determinants of M I 



8 



Therefore the proof of the theorem reduces to the calculation of numerous sub-determinants of M. These calculations 
are tedious but elementary. As a preparation, we observe that for n in {0, 1, . . . , N — 1} 



5o + 5i -5i 
-5i 5i + 52 -52 











-fly. 

~9n 9n 



90--- 9r, 



(A2) 



by repeatedly adding all other columns to the first one and pulling out factors. Repeated use of the linearity of the 
determinant in the last column of (| A2|) yields 



5o + 5i -5i 
-5i 5i + 52 -52 











~9n-l 
"An— 1 5n-l + 9n 



5o • ■ .5™^ ft *> 



(A3) 



i=0 



another identity which will be used frequently. Finally we introduce the abbreviation q = exp(— ik). 
We turn now to the calculation of the lowest order coefficients in (|Aip . 



Calculation of det M : 



detM 



fliv-i(2- q- q x ) 



3jv-i(1 - q *) 



5o + 5i -51 
-5i 5i + 52 -52 



g N -i(l - q) 

by adding all columns to the first column and subsequently all rows to the first one 



—5^-2 

5^-2 5JV-2 + 9N-1 



g N -i(2-q-q : ) 



9o + 5i -5i 
—9i 5i + 52 -52 







— flJV-2 
"fliV-2 5^-2 



by pulling out a factor and eliminating the entry in the upper right and lower left corners, 



9N-i(2-q-q 1 )go 



9i + 52 -52 
-52 52 + 53 -93 







~9N-2 
-9N-2 gN~2 



by adding all columns to the first column and subsequently all rows to the first one and pulling out a factor, 

= 5o ■ • -9N-i{2 -q- q^ 1 ) 

by applying (|A2|) . 
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Calculation of the (N — 1) x (N — 1) sub -determinants: 



Let < n < N — 1. We consider computing the sub-determinant of M where row and column n are deleted. 

-l 



det M' 



gN-i+go -go 



-gN-iq 



j n -i+g n 

1 

g n +i + g n +2 



-gN-iq 

gN-i+go -.9o 
-.90 



q gN-i 



~9n-l 

9n-l gn-1 + 9n 

gN-i+ga -.9o 

-.90 



.9^-2+3^-1 

9n+l + gn+2 —9n+2 
~9n+2 



9 N -i 





,9o + gi -,9i 
-9i 





~9N-2 

-gN-2 gN-2 + gN-i 
-.90 



93JV-1 



•• 

-gN-2 





— 9N-3 

-gN-3 gN-3 + gN-2 







' • -gN-2 





by expanding in the first and last column of the matrix. It is not hard to see that the determinants in the terms 
proportional to q and q~ 1 vanish: The corresponding matrices can be brought to a form where they contain a zero 
column by simple column operations. The remaining determinants can be treated using (|A3j) : 



N-l 



gN-igo ■ 



9n-i + 9 i 1 9n+1 ■ ■ ■ 9N - X [ 9 * 1 



i=0 



- 9n-i9o ■ ■ - 9n 
N-l 

go ■ --9N-1 X! Si 1 

i=0 



53 9i 1 9n+l 



■gN-2 



\i—n+l 
/ N-2 N 

E ^ 

\i=n+l / 



The cases where n — and n = N — 1 have to be treated separately, either by an explicit calculation or by appealing 
to the symmetry of the problem under cyclic permutations of go . . . gN-i- They yield the same result. 



Calculation of the (N — 2) x (iV — 2) sub-determinants: 

The calculation of the (N — 2) x (N — 2) sub-determinants proceeds analogously to that in the last paragraph. Let 
0<n<m<N — 1 and consider the sub-determinant of M_ where row and column m and n are deleted. Again we 
start by expanding linearly in the first and the last column: 



det M' = 

^=nm 
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9N-1 + 9o -flo 
~9o 



9n+l + 9n+2 ~9n+2 
~9n+2 



~9n-l 
~9n-l 9n-l + 9n 

9m+l + 9m+2 —gm+2 
~9n+2 



~9N-2 

-g m -i g m - 1 + g n 



q 9n-i 



ffjv-i+.9o -.9o 
~.9o 



-gjv-2 
-9n-2 gN-2 + gN-i 

l 



99 N -i 







'•■ 

-9N-2 



-go 
'■• 



' ■ —9N-2 





9n-i 



go + .9i -51 
-5i 







—9N-3 

-9ns gN-3 + gN-2 



Again it is not hard to see that the determinants in the terms proportional to q and q 1 vanish. By means of (|A3 
we get 



9n-i9o ■ ■ ■ g n 5jv-i + X] g i J 



9n-i9o ■■■gn 



J2 9* 1 



• • 5m I 2^ 9i I 5m+l • • • 9N-X I 2^ 5 * 

\i=n+l / Vi=m+1 y 

2^ 9i g m +i--.g N -2 2^ - 9 * J ■ 

\i=n+l ) \i=ra+l / 



Again, the case where the first or the last row and column get deleted have to be considered in a separate calculation 
or treated by symmetry arguments. The same result is obtained in these cases. 

When summing over all such determinants, by summing over pairs (m, n), the result further simplifies: 



V detM' 



5o • • ■ 5w 



(A4) 



0<m<n<N-l 



0<i<j<N-l 



Appendix B: Calculation of expectation values 



In the present appendix we collect some calculations regarding expectation values of the coefficients of the dispersion 
relation. 

To make the notation more simple, we will not distinguish between the random variables {1^} and a sample {k} of 
the random process anymore. The context will hopefully make the distinction clear. 
We start with the lowest order coefficient. We want to calculate 



< c2 > - = 



l 2 



(Bl) 
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It is convenient to go over to new random variables 



A„ := (B2) 



These are Gaussian distributed with expectation zero and spread one. Also introducing 

n—l n—1 



c i: =5>, c 2 = ]TA?, f=| (B3) 

n=0 n=0 

we can rewrite the expectation value as 

[t\ = ( c 2 d 2 + 2Nlgd + N 2 l 2 \ _ / c 2 S 2 + 2N Cl S + N 2 \ 

\F/ ~ \Nc 2 d 2 +2Nl Cl d + N 2 l 2 / ~ \Nc 2 5 2 + 2Nc 1 5 + N 2 / ' ( ' 

We were not able to compute this expectation value exactly. We can however expand it in a series that converges for 
small variance S. To that end, we Taylor-expand in S, obtaining 



+ (Ac 4 - <oc\c 2 N + c 2 2 N 2 ) 



S^ 4 



- 4 (2c? - 3c3 C2 7V + c lC 2 2 N 2 



iV 



+ (16c? - 28cjc 2 N + 13c 2 c 2 N 2 - c 3 N 3 ) (J?j + O (5 7 ) . 



The terms in this series are such that their expectation values can be computed through tedious but straightforward 
calculations. We use 

<(A,) fe ) = |1 9#WM n ?A°*t> (A^HOfor^j (B6) 



_^ 2 | r (fe±i) if k even 



to compute 





= TV 


(c 2 ) 


= N 


<ci> 


= 3N 2 


<c?c 2 ) 


= N(N + 2) 


(cl) 


= N(N + 2) 


(c?) 


= 15N 3 


<cfc 2 ) 


= 3N(N 2 + N + 3) 




= N(N 2 + 6N + 8) 




= N(N 2 + 6./V + 8) 



(B7) 



Monomials in Ci , c 2 that contain odd powers of the A^ vanish. Taking the expectation value of the above Taylor series 
then gives 



(-iV 6 + 7N 5 - UN 4 + 260N 3 - 2b2N 2 ) (J^j + 0{5 8 ) 



(B8) 
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A remarkable aspect of this series expansion of the expectation value is that term by term, the limit N — > oo is well 
defined, and suggests 

,. /f\ „ d 2 d A d 6 1 

lim ( = ) = 1 - — + — - — + ... = — . (B9) 

N^oo\p / P l A I 6 i + & v 1 

While the above calculations do not constitute a proof in the strict sense, we are quite confident that it is correct, 
since it is backed up both by heuristics and by numerical simulations. This is discussed in the main text. 
For this term, we can even easily compute the variance with the same methods. We first have to expand 



/ c 2 <5 2 + 2N Cl 5 + N 2 \ 1 n . 2 sA / 5 \ 



2 



^3 



+ 3(c 2 N 2 -4c 2 c 2 N + 34)^ 4 (BIO) 
- 4 (3c 2 Cl N 2 - 8c 2 clN + 5c?) (J^j 

+ (-4clN 3 + 42c 2 c 2 N 2 - 82c 2 c 4 A + 44c?) (J^j ' + O (5 7 ) 

again, the expectation value of each of the terms can be evaluated and the limit TV — > oo taken. We will only give the 
result: 

lim l(l) 2 ) = l-4 + 4-4 + ... = ^- 2 . (Bll) 

Therefore we obtain the very simple result that 

lim (c 4 ) - lim (c 2 ) 2 = 0. (B12) 

The situation for the coefficient t 2 of \k\ 4 is more murky. Again we change variables and obtain, employing the same 
notation as above, 

[ ' i<j y ' i<j (B13) 

1 I 2 ( Cl< 5 + A^) 4 
~ 12Nc 2 S 2 + 2 Cl S + N' 

Again, we try an expansion in 5. We use 

^ + N ? -N-2 Cl S + 3 ^; T C2N h 2 + Q(^ (B14) 



(2ci<5 + c 2 S 2 + N) 3 ~ ' JV 

and 



(ciJ + JV) 



4 



= n 3 + 2 Cl n 2 5 + (2cjn - c 2 n 2 ) S 2 (B15) 



c 2 5 2 + 2 Cl 5 + N 
to find after some calculation 

^ = _ S + l^ -4 ^ + N ~ 6)62 + ° {S3) (B16) 

The problem with this expansion is that, at least term by term, the N — > oo limit can not be taken: The term 
proportional to S 2 diverges linearly with N. This is despite some cancellation between the two terms that comprise 
C. We will discuss the implications of this in the main text. 
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Appendix C: Discrete time, continuous space 



As metioned in the introduction, we can also apply our calculations to the situation in which time is discrete and 
space is continuous and homogeneous, by exchanging the role of space and time in all formulas. The action in this 
case would read 



S= f dxJ2[9n(d + ^ n -^' n ) 2 ] 

■1 ra6Z 



with the prime denoting the spatial derivative of the field 4> n (x) . Determination of the equations of motion proceeds 
in parallel with the main text. The wave-ansatz is now 

N-l 

4>n >{ x ) = c n expi(Lzu! — kx), L — l„. (CI) 
The calculation of the main text then gives the dispersion relation 

fc 2 H = c 2 uj 2 + £ 2 lj 4 + . . . , (C2) 

with c 2 and £ 2 exactly as in ([15]), and (|16|) . We can invert this series to get to the more familiar form 

w 2 (fc) = lfc 2 - Jfc 4 + .... (C3) 

We have not analyzed the expectation value (c~ 2 ) in detail, but one can argue that since the distribution of c 2 is 
sharp in the limit N — > oo one should have 

jrj^>=tej<r'=( i+ £) 2 - (C4) 

i.e. we obtain a positive correction to the continuum phase velocity. Also, the coefficient of fc 4 comes with a negative 
sign, thus presumably leading to an increase in group velocity, as well. 



L. Bombelli, J. H. Lee, D. Meyer and R. Sorkin, "Space-Time As A Causal Set," Phys. Rev. Lett. 59 (1987) 521. 
R. Gambini and J. Pullin, "Nonstandard optics from quantum spacetime," Phys. Rev. D 59 (1999) 124021 
[arXiv:gr-qc798 09038 . 

H. Sahlmann and T. Thiemann, "Towards the QFT on curved spacetime limit of QGR. I: A general scheme," Class. Quant. 
Grav. 23 (2006) 867 [arXiv:gr-qc/0207030l . 

H. Sahlmann and T. Thiemann, "Towards the QFT on curved spacetime limit of QGR. II: A concrete implementation," 
Class. Quant. Grav. 23 (2006) 909 [arXiv: gr-qc/0207031 1 . 

S. Bernadotte and F. R. Klinkhamer, "Bounds on length scales of classical spacetime foam models," Phys. Rev. D 75 
(2007) 024028 [arXiv:hep-ph/0610216| . 

H. Sahlmann, Coupling matter to loop quantum gravity, PhD Thesis, Potsdam University (2002) [URN: 
urn:nbn:de:kobv:517-0000602, URL : |http : //opus.kobv. de/ubp /volltexte /2005 /55/] 

N. H. Christ, R. Friedberg and T. D. Lee, "Random Lattice Field Theory: General Formulation," Nucl. Phys. B 202 
(1982) 89. 

N. H. Christ, R. Friedberg and T. D. Lee, "Gauge Theory On A Random Lattice," Nucl. Phys. B 210 (1982) 310. 
F. R. Klinkhamer and C. Rupp, "Spacetime foam, CPT anomaly, and photon propagation," Phys. Rev. D 70 (2004) 
045020 [arXiv:hep-th /0312032| . 

F. R. Klinkhamer and C. Rupp, "Photon-propagation model with random background field: Length scales and Cherenkov 
limits," Phys. Rev. D 72 (2005) 017901 arXiv:hep-ph/0506071| . 

L. Bombelli, A. Corichi and O. Winkler, "Semiclassical quantum gravity: Statistics of combinatorial Riemannian geome- 
tries," Annalen Phys. 14 (2005) 499 [arXiv: gr-qc /0409006] . 
C. Kittel, Introduction to Solid State Physics, Wiley (1995) 
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